Evidence for non-Gaussianity in the CMB 
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ABSTRACT In a recent Letter we have shown how COBE-DMR maps may be used to disprove 
Gaussianity at a high confidence level. In this report we digress on a few issues closely related 
to this Letter. We present the general formalism for surveying non-Gaussianity employed. We 
present a few more tests for systematics. We wonder about the theoretical implications of our 
result. 
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1 Introduction 

It is hard to overemphasise the importance of Gaussianity in theories of structure 
formation. Under the assumption of Gaussianity calculations become much simpler, 
and only the power spectrum of fluctuations has to be computed. Furthermore it is 
thought that Gaussianity is a hallmark of the inflationary paradigm Q. A simple 
argument for this hinges on the fact that small perturbations during inflation sat- 
isfy an harmonic oscillator equation, at the relevant scales for structure formation. 
The ground state of an harmonic oscillator has a Gaussian wave function. Hence 
although an unperturbed background is the classical ground state, quantum the- 
ory forces the existence of "zero-level" fluctuations with Gaussian statistics. These 
quantum fluctuation are the seeds of structure, according to inflation, and one iden- 
tifies quantum probability with the classical statistics of these macroscopic seeds. 
This simple argument linking inflation and Gaussian statistics is bypassed in non 
minimal models of inflation (|]. 

In a recent Letter [|| we showed how Gaussianity could be disproved at a high 
confidence level, using COBE-DMR maps. Here we review that Letter and report 
on work in progress expanding || . In Section 2 we describe a generalisation of the 
method employed in which incorporates inter- £ correlators. In Section 3 we 
review how the method was applied to the data in j|. In Section 4 we review the 
checks for systematics considered in || and present a few extra tests. 

We conclude with a few remarks on recent work ||, ||, |?j which comments on 
our result. 
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2 The statistics 



The statistics used in |j| are a subset of a general class of statistics, to be presented 
in [Ej. They are inspired by || [nj. The idea is to extract from a laboriously 
produced CMB map all the relevant information, with as little theoretical prejudice 
as possible. If we believe inflation is the answer, then one may simply reduce a 
map to a fit for a set of cosmological parameters. A less omniscient approach is 
to assume only Gaussianity, and concentrate on an unbiased estimate of the CMB 
power spectrum (see fl2|| for an example). The ultimate open mind would not 
assume Gaussianity, but try to extract from the map the whole set of correlation 
functions characterising the most general random process. 

It is fair to keep one prejudice: statistical isotropy. The idea is therefore to ex- 
tract from a map with N pixels ./V — 3 rotationally invariant independent quantities. 
We work in the spherical harmonic representation a m . To construct a an ri-linear 
invariant one takes the tensor product of n ATgs 

(AT tl (g> AT e2 (g> ■ ■ ■ (g> AT e J(n) = 

nil 7712 f^n 

x Y timi (n) <g> Y e2Jn . 2 (n) ® • • • ® Y tn7nn (n) (1) 

One is interested in rotationally invariant quantities. These can be trivially obtained 
if one rewrites the tensor product in terms of the total angular momentum basis. 
The coefficient of the singlet will be the higher order invariant we are looking for. 

To illustrate the technique we shall work out a few examples. Let us first 
construct the rotationally invariant bilinear. For this we take the tensor product 
AT £l (g> AT l2 : 

(AT tl <g> AT fe )(n) = ^2 a limi a l2m2 Y limi (n) ®Y/, 2m2 (n) (2) 

mim2 

One can now use the angular momentum addition formulas to find the the coefficient 
of the singlet. Defining 

Y limi (n) ® Y kma (n) = ^){LM\U mi m 2 )Y LM {n) (3) 

LM 

we want the L = term. This implies that t\ — £2 — £ and from the condition 
M = we have mi + mi = 0. So the coefficient of the singlet is 

(—l) e ~ m 

m V i 

Up to normalisation this is the simplest quadratic estimator of the power spectrum. 

The next simplest case is the cubic invariant. Let us first restrict ourselves to 
one "ring" , i.e. fixed £: 

(ATi ® ATg <g> AT e )(n) 
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= ^ dinn aim 2 fl|m 3 ^imi (n) ® Y(? 2m2 (n) ® Y^ m3 (n) 

mim2m3 

= ^ atm 1 «<m 2 fl<m 3 ^ (LM\££m 1 m 2 ) Y LM (n) ® Fi m3 (n) 

mim2inj LM 

a emi a em2 a em3 (LM\£em 1 rn2)(L'M , \L£Mm 3 )Y L ,M'(n) 

mxmimz LL' MM' 

Once again, we want the singlet. From L 1 = M' = we get L = I and M = — 777.3 • 
This leaves us with 

a£m 1 atm 2 a em 3 (i - m 3 \££mim 2 ) (00\££ - m 3 m 3 ) 

ra\ra 2 ra 3 

1 , ^ / * £ £ \ 

£+1) \ m l m 2 m 3 J 



From the symmetry properties of the Wigner 3 J coefficients we immediately see 
that for I odd this quantity is identically zero. 

We are also interested in relating power between is. Consider then (ATf_i <g) 
ATe ® ATf + i)(n). From the same manipulations one finds 

J e = Y a t-im 1 ai+i m2 ai> m3 (i ~ m 3 \£ + 1£ - lmxm 2 ) (00\£i - m 3 m 3 ) 
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1 e ^ ( i-1 £+1 £ x 

z \ / . a £rni a lm2 a Zrn,3 I I %i+m2+m3,0 

1 ^— » \ 777 1 777, o 777, t / 



J (21+ IV ' ^ 1 A V m i TO 2 m 3 

V \ / mi 7719 7713 



(5) 



This procedure can be used to find all the invariants at any order, within a given 
multipole (I 1 ) and relating different multipoles (J e ). These may then be divided by 
the appropriate powers of the Cg in order to make them dimcnsionlcss, and suitably 
normalised, as was done with If in ||. 

This method should produce the full set of independent invariant quantities in a 
set of a l m . There should be 2i—2 such quantities, for each £, plus 3 inter-£ invariants, 
for each pair of is. The power spectrum measures how much power there is on a 
given scale £. The I describe how the power is divided between the various m 
modes, for a given i. This encodes preferences for shapes within a given multipole. 
The J e measure correlation between the orientations of preferred shapes in adjacent 
multipoles. 



3 DMR bispectrum 

We now summarise the application of the simplest of these statistics to DMR. 
We will be testing the inverse noise variance weighted, average maps of the 53A, 
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Fig.l The vertical thick dashed line represents the value of the observed if. The solid line is 
the probability distribution function of if for a Gaussian sky with extended galactic cut and 
DMR noise. 

53B, 90A and 90B COBE-DMK channels, with monopole and dipole removed, at 
resolution 6, in ecliptic pixelization. We use the extended galactic cut of [jl3|, and 
[ [Til to remove most of the emission from the plane of the Galaxy. We apply our 
statistics to the DMR maps before and after correction for the plausible diffuse 
foreground emission outside the galactic plane as described in jl5| and |Tq| . To 
estimate the ifs we set the value of the pixels within the galactic cut to and the 
average temperature of the cut map to zero. We then integrate the map multiplied 
with spherical harmonics to obtain the estimates of the ae m s and compute the If 
from these. 

We have used Monte Carlo simulations to find the distribution of the estimators 
if as applied to Gaussian maps subject to DMR noise and galactic cut (see Fig.l). 
These distributions are very non-Gaussian. In principle this would complete the 
theoretical work required for converting the observed if (which we also plot in Fig.l) 
into a statistical statement on Gaussianity, but we proceed further by defining a 
new "goodness of fit" statistic as follows. 

We constructed a tool similar to the % 2 (often used for comparing predicted and 
observed Ci spectra) but adapted to the non-Gaussian distributions P{lf): 

X " = ^E X * = ^£(-21og^(J?) +&), (6) 
e e 
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where the constants /3e are dehncd so that for each term of the sum (Xf) = 1. The 
definition reduces to the usual X 2 for Gaussian Pi. 

This is a suitable definition. As an illustration let us first approximate the 
distributions of the If by P(lf) = 2(1 — if) — a good approximation for £ around 
10. Then X 2 = -21og(l - if). Like the standard X 2 one has < X 2 < 1 for 
observations close to the peak of the distribution, here at if — 0. Indeed X 2 (0) = 0. 
However the peak of P(If) is far from its average, and so the standard X 2 would 
produce X 2 — at the wrong observation. For observations far from the peak of the 
distribution (but subject to the constraint if < 1) X 2 goes to infinity. In contrast 
the standard X 2 would always remain finite. 

The proposed X 2 therefore does for these non-Gaussian distributions what the 
usual X 2 does for normal distributions. We build a X 2 for the COBE-DMR data 
by means of Monte Carlo simulations. We proceed as follows. First we compute the 
distributions P(lf), for I = 2, 18, for a Gaussian process as measured subject to 
our galactic cut, and pixel noises. These P(lf) were inferred from 25000 realizations 
(see Fig.l). From these distributions we then build the X 2 as defined above, taking 
special care with the numerical evaluation of the constants fy. We call this function 
-X-COBE- We then find its distribution F{X^ OBE ) from 10000 random realizations. 
This is very well approximated by a \ 2 distribution with 12 degrees of freedom. If 
all P(If) were as in the analytical fit above, we could conclude that we successfully 
measured an effective number of useful invariants equal to 6. This is less than the 
number of invariants we actually measured (10) and this is simply due to anisotropic 
noise and galactic cut. However, had we used a standard x 2 statistics the effective 
number of useful invariants would be only 3. 

We then compute X^ OBE with the actual observations and find X^, OBE = 1.81. 
One can compute P(X^ OBE < 1.7) = 0.98. Hence, it would appear that we can 
reject Gaussianity at the 98% confidence level. 

4 Systematics 

Given the nature of this result checking for systematics has been the central aspect 
of our work. We checked a large number of effects, related to foreground emissions, 
pixelization effects, spurious offsets induced by the galactic cut, the cut itself, the 
underlying power spectrum, etc, etc. These checks are reported in ||, 

It is important to stress that the confidence level quoted above is a lower limit. 
It corresponds to the worse case obtained, within the systematic space surveyed. It 
assumes a conspiracy theory with all systematics lined up so as to take the blame for 
the observed non-Gaussianity. If one does not consider this worst case scenario, but 
takes other data sets differently treated for systematics, in most cases we obtained 
a confidence level for rejection in excess of 99.5%. 

Checking for systematics is always open ended, and in these Proceedings we 
merely present one new check we have recently performed. It concerns the so-called 
systematic templates. The procedure leading from the time series to a CMB map 



5 




4 8 12 16 3D 4 8 1 2 1 6 20 



Fig. 2 - The effect of adding and subtracting templates of systematics due the instrument 
magnetic susceptibility, in directions px, Pr, and px> an d the effect of potential unknown 
effects at spin frequency. 

produces, as well as the map, an estimate of the systematic effects possibly plaguing 
the final product. These "systematic templates" for DMR are well documented 
(fTi]]), and display strongly non-Gaussian structures, tracing the DMR scanning 
patterns. 

The systematic templates have negligible power. The worst effect in the worst 
channel has a rms of about 6 \iK at the 95% confidence level. It is unlikely that these 
effects could corrupt power spectrum estimates. Nonetheless it is well known that a 
non-Gaussian pattern with negligible power may visually stand out over a Gaussian 
map with much larger power. Similarly it could happen that these systematics, 
while irrelevant for the purpose of power spectrum estimation, could be responsible 
for the observed nonGaussian bispectra, derived from DMR maps. 

In order to address this problem, we subjected the systematic templates to two 
tests. Firstly we computed the if spectra for the templates. The resulting If 
are well outside the Gaussian prediction, but they do not correlate with the DMR 
observed If. Secondly, we added or subtracted these templates enhanced by a factor 
of up to 4 to DMR maps. The effect on the If spectrum was always found to be 
negligible. This shows that the systematic effects documented in have not only 
negligible power, but also negligible effect on the higher order statistics which we 
have studied. 

To be more specific we have applied the above tests to templates for systematics 
in 53A, 53£>, 90A, and 90B, separately. This is the sensible thing to do, given that 
the templates are highly correlated from pixel to pixel We have considered the 
effect of instrument susceptibility to the Earth magnetic field; any unknown effects 
at the spacecraft spin period; errors in the calibration associated with long-term 
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drifts, and calibration errors at the orbit and spin frequency; errors due to incorrect 
removal of the COBE Doppler and Earth Doppler signals; errors in correcting for 
emissions from the Earth, and eclipse effects; artifacts due to uncertainty in the 
correction for the correlation created by the low-pass filter on the lock-in amplifiers 
(LIA) on each radiometer; errors due to emissions from the moon, and the planets. 
Presumably what is usually mentioned as De-striping goes under the removal of 
artifacts in the calibration. DMR did not seem to have a serious striping problem, 
but the problem was addressed none the less. 

In Fig. 2 we plot the result of adding templates associated with effect of the 
magnetic field of the Earth. The instruments' magnetic susceptibility, and the 
emissions from the Earth have by far the strongest effects. The effect on the if 
spectrum is always negligible, but when present occurs at scales around I = 6. 

5 Cautionary remarks and a digression 

By now two other groups have reported results similar to ours, albeit making use of 
different methods ||. According to skeptics, this may merely reflect a change in 
the psychological prior, triggered by our work. More seriously one should remember 
that the work performed by us and by these groups makes use of the same data set. 
Therefore this work provides an independent confirmation of our analysis of the 
DMR maps, but not an independent confirmation of the result itself. In particular 
we feel that the issue of systematics, and foreground contamination, will only be 
clarified further when an independent all sky data set becomes available. 

If indeed our result is due to cosmic emission, then a number of fascinating the- 
oretical issues are raised (see |?]]). Clearly the minimal inflationary models cannot 
be right. On the other hand it is not obvious that the main competitor to infla- 
tion, topological defects, could explain this type of non-Gaussianity. Topological 
defects are non-Gaussian, but in ways which are often more subtle than commonly 
thought. Computing with defects is prohibitively expensive, and predicting a set of 
If distributions in defect models is well beyond current computer technology. 

An interesting possibility was recently proposed by Peebles Q . This is an isocur- 
vature model in which the underlying fluctuations are not a Gaussian random field, 
but the square of a Gaussian random field. The model is based on non minimal 
inflation, but produces fluctuations radically different from minimal inflationary 
fluctuations. The main advantage of this model over defects is that, while not triv- 
ial, it is easy enough to compute with it. In particular it is feasible, and topical, 
to repeat the exercise we have performed in || using Peebles theory. This model 
could well produce a better fit to the DMR bispectrum than Gaussian theories. 
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